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ABSTRACT: 

The  Flood  Event  Assessment  Tool  (FEAT)  is  a  prototype  geospatial  modeling  tool  that  uses  river  and 
stream  gage  data,  landscape  digital  elevation  models  (DEMs).  main  and  secondary  channel  centerlines  or 
cross  sections  to  generate  a  geospatial-based  flood  surface.  The  primary'  objective  in  developing  FEAT 
was  to  improve  the  accuracy  and  decrease  the  time  and  effort  required  to  determine  dynamic  flood  sur¬ 
faces  and  provide  information  for  decision-making.  This  guide  provides  an  overview  of  FEAT,  informa¬ 
tion  about  the  algorithms,  software  and  hardware  requirements,  installation  procedures  and  operation, 
user  input  requirements,  and  examples  of  model  output  products. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes. 
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All  product  names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not 
to  be  construed  as  an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 
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District. 

The  report  was  prepared  by  Mr.  Jerrell  R.  Ballard,  Jr.,  and  Dr.  Margaret  Rose 
Kress  of  the  Environmental  Characterization  Branch  (ECB),  Natural  Resources 
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Introduction 


Existing  methods  for  estimating  the  extent  of  flooding  in  the  Lower 
Mississippi  River  Valley  primarily  consist  of  interpretation  of  satellite  data. 

These  methods  are  labor-intensive  and  often  are  restricted  due  to  extensive  cloud 
cover  during  a  flood  event.  Another  limitation  includes  the  inability  to  “project” 
flood  extent  for  future  events  or  Corps  project  developments. 

The  Flood  Event  Assessment  Tool  (FEAT)  is  a  prototype  geospatial 
modeling  tool  that  uses  river  and  stream  gage  data,  landscape  digital  elevation 
models  (DEMs),  and  main  and  secondary  channel  centerlines  or  cross-sections  to 
generate  a  geospatial-based  flood  surface.  The  primary  objective  for  the 
development  of  FEAT  was  to  improve  the  accuracy  and  decrease  the  time  and 
effort  to  determine  dynamic  flood  surfaces  and  provide  information  for  decision 
making.  The  conceptual  algorithm  for  the  geospatial  flooding  was  initially 
developed  for  the  Flood  Impact  Support  Tool1  (FIST)  system  and  was 
additionally  documented  by  Ballard  and  Kress.2 

This  guide  provides  an  overview  of  FEAT,  information  about  the  algorithms, 
software  and  hardware  requirements;  installation  procedures  and  operation;  user 
input  requirements,  and  examples  of  model  output  products. 

This  modeling  procedure  was  developed  by  the  U.S.  Army  Engineer 
Research  and  Development  Center,  Waterways  Experiment  Station,  Vicksburg, 
MS,  for  the  Hydrography  and  Hydrology  Division,  Vicksburg  District,  U.S. 
Army  Corps  of  Engineers. 


1  Ballard,  J.  R.,  Jr.,  and  Kress, M.  R.  (1999).  Flood  Impact  Support  (FIST)  User's  Manual 
and  Technical  Documentation.  Instruction  Report  EL-99-1,  U.S.  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  MS. 

2  Ballard,  J.  R.,  Jr.,  and  Kress,  M.  R.  (1998).  “Cost-Distance  Flood  Surface  Generation 
for  Economic  Impact  Assessment.”  Proceedings  for  the  1 8th  Annual  ESRI  User 
Conference,  San  Diego,  California,  July  27-31,  1998. 
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2  Components  of  FEAT 


FEAT  combines  spatial  analysis  techniques  and  high-resolution  geospatial 
data  to  determine  land  area  inundated  under  specified  flood  conditions,  FEAT 
has  four  components: 

•  A  fully  functional  commercial  off-the-shelf  (COTS)  geographic 
information  system  (GIS)  (ArcView®3.1.1). 

•  A  high-resolution  geospatial  database  for  the  area. 

•  The  flood  impact  model. 

•  A  custom-designed  user  interface. 

The  GIS  serves  as  the  database  manager,  provides  access  to  a  complete  range 
of  spatial  analysis  capabilities,  and  establishes  the  graphic  display  environment. 
ArcView®firom  Environmental  Systems  Research  Institute  (ESRI)  is  the  GIS 
used  for  FEAT. 

The  high-resolution  databases  depict  quantitative  information  on  the  physical 
terrain  surface,  water  elevations  at  stream  recording  gages,  and  land-cover 
(vegetation  types)  conditions  in  selected  study  areas.  Required  database 
components  are  documented  in  Chapter  4,  “Database  Requirements.” 

The  flood  impact  model  has  been  designed  and  developed  to  run  in 
conjunction  with  and  utilize  the  capabilities  of  ArcView®.  It  includes  programs 
written  in  Avenue,  the  programming  language  within  ArcView®.  These  programs 
are  listed  in  Appendix  A. 

The  user  interface  is  written  using  the  Avenue  language.  The  interface  helps 
the  user  select  appropriate  data  types  for  the  model  and  establishes  user  set 
parameters. 
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3  Hardware,  Software,  and 
Training  Requirements 


FEAT  requires  specialized  computer  hardware  and  software  for  proper 
system  operation.  Simulating  flood  events  is  computationally  intensive  and 
requires  a  high  performance  personal  computer  platform  for  optimum 
performance.  To  achieve  this  level  of  performance,  the  minimum  hardware 
requirements  are  a  computer  running  Windows  98  with  at  least  a  400-Mhz 
Pentium  II  CPU,  CD-ROM  drive,  6  gigabytes  of  hard  disk  space,  and  a  16-bit, 
17-in.  color  graphics  monitor. 

The  GIS  used  in  FEAT  is  ESRI  ArcView®,  with  the  Spatial  Analyst 
Extension.  ESRI  Avenue  code  was  utilized  for  the  algorithms  in  FEAT.  Because 
of  this,  the  minimum  software  requirements  are  ESRI  ArcView®  Version  3.1.1, 
and  ESRI  Spatial  Analyst  Extension  Version  1.1. 

It  is  recommended  that  the  end  user  of  FEAT  be  trained  and  well  acquainted 
with  the  operation  of  ArcView®  and  have  a  technical  understanding  of  geographic 
information  systems  concepts  (GIS)  and  hydrology.  ESRI  provides  several 
training  courses  during  each  year  at  locations  nationwide. 
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4  Database  Requirements 


Specific  data  types  are  required  by  FEAT.  Each  type  of  data  are  briefly 
described  below. 

FEAT  is  designed  to  use  data  contained  in  the  GIS  database  and  other 
usersupplied  parameters.  All  geospatial  data  used  in  FEAT  must  be  projected  to 
the  Universal  Transverse  Mercator  (UTM)  projection,  North  American  Datum 
(NAD83).  The  required  GIS  data  layers  are  as  follows: 


Table  1 

Required  GIS  Data  for  FEAT 

Name 

Description 

Elevation 

Basin  ground  surface  elevation  model  (Required) 

Water  Gages 

Location  of  gages,  and  attributes  of  water  surface  elevations  (Required) 

Cross  Sections 

Channel  cross  sections  and  attributes  of  water  surface  elevations 
(Optional) 

River  Center  Lines 

River/Channel  center  lines,  attributed  with  a  from  and  to  gage  (Required) 

Secondary 

Channels 

Secondary  channel  center  lines  (Optional) 

Landcover 

Landcover  (Optional) 

Elevation 

The  “elevation”  file  is  a  grid  file  containing  the  digital  elevation  model 
(DEM)  of  the  basin.  This  DEM  must  be  able  to  accurately  represent  the 
topography  of  the  basin.  This  implies  that  all  channels,  streams,  and  levees  must 
be  represented  with  an  adequate  grid  cell  resolution. 

The  DEM  can  be  processed  from  digitized  from  1:24,000  U.S.  Geological 
Survey  (USGS)  topographic  maps  and  supplemental  field-surveyed  elevations. 
The  DEM  can  also  be  obtained  from  the  USGS  as  DEM  level  II  products.  As  a 
general  rule,  for  any  project,  the  cell  resolution  must  be  at  least  four  times 
smaller  than  the  width  of  the  smallest  river  channel  in  the  data  set.  If  the  grid  cell 
resolution  is  too  large,  small  channel  geometry  may  not  be  represented 
sufficiently  in  the  DEM  and  thereby  result  in  inaccurate  modeling  of  the  flooded 
area,  especially  backwater  flooding.  The  DEM  vertical  units  must  be  represented 
as  meters  above  mean  sea  level. 
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Water  Gages 

The  “water  gages”  file  is  a  point  feature  file  containing  the  geographic 
location  and  the  unique  identification  number  of  each  river  gage.  Also,  at  each 
river  gage  there  must  be  at  least  one  water  surface  elevation  associated  with  each 
river  gage  location.  Table  2  provides  an  example  of  possible  water  gages  point 
file  attributes  where  the  water  surface  elevations  are  provided  in  feet  above  mean 
sea  level.  These  water  surface  elevations  can  be  obtained  from  either  direct 
readings  or  model  outputs. 


Table  2 

Water  Gage  Required  Items 

Example) 

Gage  ID 

Name 

Two-Year  (ft) 

Ten Year (ft) 

1 

Anguilla 

98.0 

100.2 

2 

Holly  Bluff 

95.3 

98.2 

Cross  Sections 

The  “cross  sections”  file  is  a  line  feature  file  containing  the  channel  cross 
sections  with  attached  water  surface  elevation  attributes.  These  cross  sections  can 
be  obtained  typically  from  Hydrologic  Engineering  Center  (HEC)  model  outputs. 
This  file  can  be  used  optionally  with  FEAT. 


River  Center  Lines 

The  “river  center  lines”  file  is  a  line  feature  file  containing  river  and  channel 
center  lines.  Each  line  is  attributed  with  an  upstream  gage  identification  number 
and  a  downstream  gage  identification  number.  These  attributes  allow  FEAT  to 
identify  the  proper  up-  and  downstream  water  elevations  at  run-time.  Table  3 
provides  an  example  of  the  attributes  for  one  center  line.  In  this  example  the  river 
center  line  runs  from  gage  1  (Anguilla)  to  gage  2  (Holly  Bluff)  according  to 
Table  2. 


Table  3 

River  Center  Line  Required  Items  (Example) 

Center-line_ID 

From 

To 

1004 

1 

2 

Secondary  Channels 

The  “secondary  channels”  file  is  a  line  feature  file  that  represents  important 
secondary  channels  in  the  flood  event.  These  are  unattributed  line  segments  that 
are  directly  attached  to  a  river  center-line  vertex.  During  execution,  FEAT 
determines  the  water  surface  elevation  at  the  center-line  vertex  and  associates  it 
with  the  entire  secondary  channel.  A  schematic  example  of  this  is  shown  in 
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Figure  1.  The  secondary  channel  file  is  useful  in  obtaining  proper  backwater 
flooding  in  small  areas  with  inadequate  delineated  channels. 


Figure  1 .  Main  channel  center  line  and  secondary  channel  example 


Landcover 

The  “landcover”  file  is  a  grid  file  describing  cleared,  wooded,  and  urban 
areas  in  the  basin.  This  grid  file  is  used  after  a  flood  surface  is  generated  to 
determine  which  types  of  landcover  are  impacted  by  the  flood  event. 
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5  FEAT  Operation 


The  FEAT  extension  in  Arc  View®  contains  a  simple  graphic  interface  that 
requires  inputs  from  the  user.  The  following  text  describes  the  operation  of  the 
extension,  interface,  user  inputs,  and  examples  of  its  operation. 


Activating  the  FEAT  extension 


Load  the  extension  using  the  extensions  dialog  under  the  File  pull-down 
menu.  Select  and  load  the  FEAT  extension.  The  user  should  see  a  dialog  box  that 
looks  similar  to  Figure  2.  Find  where  the  extension  is  listed  and  place  a  check  in 
the  box.  Then  press  the  OK  button  on  the  Extensions  dialog  to  load  the 
extension.  The  extension  will  load  itself  and  the  Spatial  Analyst  extension  into 
Arc  View  and  is  now  ready  for  use. 


Be  Project  Window  Help 


raiaSaiIliyiS2l! 


New 


Available  Extensions: 


!®>4 


T  ables 

IE 


_|  Legend  T  ool  jj|] 

_J  MrSID  Image  Support 
_J  NITF  Image  Support 

Project  Flood  Impact  Assessment 
_|  Report  Writer 
_J  Script  Editor  Utilities 
_|  Spatial  Analyst  _ O 


OK 


*1 


J 


Cancel  | 


_ Reset 

Or  Make  Default 


Charts  About: 


Layouts 


Scripts  zJT 


Figure  2.  Enabling  the  FEAT  extension  within  ArcView 
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Starting  the  Flood  Event  Simulation 


Once  the  required  data  layers  are  installed  into  a  view,  (the  ArcView®  users 
manual  explains  how  to  do  this),  the  display  should  look  similar  to  Figure  3.  The 
following  steps  explain  and  user  inputs  and  actions. 


Figure  3.  Initial  setup  of  required  data  layers 


STEP  1:  The  pull  down  menu  is  titled  ‘Flood  Impact  Assessments’  and 
select  the  one  ‘Determine  Flood  Extent’  (Figure  4).  The  extension  will  prompt 
the  user  to  indicate  the  elevation  grid  theme  to  be  used  for  calculations 
(Figure  5),  If  no  grid  themes  are  available,  the  simulation  will  terminate.  There 
will  be  only  one  elevation  grid  theme  used  during  the  simulation. 
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Figure  4.  Selecting  the  option  in  the  pull-down  menu 


Figure  5.  Elevation  grid  theme  prompt 
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Figure  6.  Prompt  for  the  Water  Gage  theme 


STEP  3:  Once  selected,  an  examination  of  the  attributes  (or  fields)  attached 
to  the  points  will  be  made  to  indicate  the  water  surface  elevations. 


STEP  4:  The  last  prompt  associated  with  the  gages,  is  to  check  the  units  of 
the  water  surface  elevations  (Figure  7).  If  the  units  are  not  in  feet,  then  the  units 
are  assumed  to  be  in  meters. 


Figure  7.  Vertical  units  prompt  for  the  water  surface  gages 


STEP  5:  The  next  step  is  to  select  a  centerlines  theme.  All  themes  are 
examined  in  the  view,  and  only  the  line  themes  are  presented.  It  is  the  user’s 
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Figure  8.  River  Center  Lines  prompt 


STEP  6:  This  step  examines  the  attributes  for  the  selected  point  theme  and 
prompts  for  item  names  where  the  beginning  and  ending  water  gage 
identification  numbers  are  stored  (Figures  9  and  10). 


Figure  10.  Center  Line  Water  Gage  Downstream  selection 


STEP  7:  To  enhance  the  water  surface  interpolation,  additional  vertices  may 
be  added  to  the  centerlines  by  increasing  the  density  of  all  the  lines.  A  prompt 
asks  for  the  density,  but  if  a  density  increase  is  not  required  then  the  density 
should  be  set  to  a  number  greater  than  the  average  density  of  the  centerline’s 
vertices  (Figure  1 1). 


Figure  1 1 .  Vertex  density  selection  for  the  stream  center  lines 


STEP  8:  If  a  center-line  theme  is  selected  for  use  in  the  water  surface 
interpolation,  the  next  prompt  determines  if  any  secondary  channels  are  to  be 
used  (Figure  12).  There  are  no  attributes  associated  with  this  line  theme,  and  the 
simulation  will  determine  which  end  point  of  the  line  connects  to  a  center  line 
and  will  assigns  a  water  surface  elevation  to  the  entire  secondary  channel. 
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Figure  12.  Secondary  channel  selection 


STEP  9:  The  next  prompt  is  for  river  cross  sections.  This  is  similar  to  the 
river  center  lines  selection.  Once  the  line  theme  has  been  selected,  an 
examination  of  the  attributes  is  made  and  prompts  the  user  again  for  the  attribute 
(item)  that  contains  the  water  surface  elevations. 

STEP  10:  In  step  10,  the  user  is  prompted  for  the  number  of  steps  to  use  in 
the  water  flow  calculations.  The  default  is  25  steps  (Figure  13).  This  may  be 
incremented  to  slightly  increase  the  accuracy  of  the  water  flow  calculations,  but 
dramatically  increases  the  run  time  of  the  simulation.  An  explanation  of  how  this 
parameter  works  and  how  it  affects  accuracy  of  the  water  flow  calculations  is 
provided  in  Chapter  7. 


Figure  1 3.  Sampling  steps  for  the  water  flow  calculations 


After  this  last  prompt,  modeling  of  the  flood  event  is  initiated. 
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Operational  Guidance 


During  testing  and  developing  of  data  sets  for  testing  the  simulation  model, 
several  operational  issues  were  noted.  The  following  is  a  general  list  of 
operational  guidance  hints: 

a.  All  river  gage  locations  (x,y  coordinates)  must  be  located  in  the  center  of 
the  river  channels. 

b.  River  centerlines  must  be  located  in  the  center  of  the  river  channels. 

c.  All  gage  readings  must  have  a  valid  water  elevation  (no  zeros  or  missing 
values)  and  set  to  an  elevation  above  mean  sea  level. 

d.  Increasing  sampling  size/rate  of  the  water  flow  calculations  (25-50)  will 
provide  better  results.  Increasing  this  rate  will  roughly  double  simulation 
times. 

e.  Simulation  times  will  increase  with  the  flooded  area  extent. 

f.  Increasing  the  density  of  vertices  of  river  center  lines  will  improve 
results.  Increasing  the  density  smaller  than  twice  the  elevation  cell  size 
will  not  improve  results. 

g.  Adding  secondary  channels  in  small  backwater  flood  areas  will  improve 
results  and  minimize  external  water  channel  influences. 

h.  In  ArcView®  set  the  project  properties  working  directory  to  a  directory 
other  than  /windows/temp.  All  grids  and  final  flood  extent  surfaces  will 
be  saved  to  the  working  directory  specified. 
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6  Diagnostic  Mode 


Two  versions  of  the  FEAT  extension  are  available  in  ArcView®.  Both  of 
these  extensions  provide  the  same  end  results,  but  one  extension  has  been 
modified  to  provide  a  diagnostic  mode  and  saves  all  intermediate  geospatial  data 
generated  during  the  flood  event  simulation.  These  intermediate  data  are  useful  in 
diagnosing  and  verifying  simulation  input  data  and  understanding  how  the 
iterative  simulation  determines  the  flood  extent  within  the  basin  or  region. 

In  addition  to  the  final  flooded  surface,  in  the  diagnostic  version  of  FEAT, 
there  are  also  seven  other  intermediate  diagnostic  grid  and  point  themes  that  can 
be  used  for  understanding  how  the  simulation  is  determining  which  land  areas  are 
flooded  and  how  water  is  routed  over  the  land. 

The  first  data  theme  produced  in  the  diagnostic  mode  is  a  grid  theme  named 
TnitSurf  (Figure  16).  This  grid  theme  is  a  grid  of  water  surface  elevations 
generated  by  using  all  the  gage,  centerline,  and  secondary  channel  water  surface 
elevations. 

The  second  data  theme  produced  in  the  diagnostic  mode  is  a  temporary  grid 
theme  with  the  name  ‘Gridxx,’  where  ‘xx’  is  a  number  generated  during  the 
simulation  (Figure  17).  This  grid  is  the  result  of  the  first  flood  event 
approximation  where  the  simulated  water  surface  elevation  is  compared  to  the 
digital  elevation  surface.  The  resulting  grid  only  shows  the  extent  of  the  flooded 
area.  Only  the  grid  cells  that  could  be  flooded  under  these  conditions  are 
indicated.  Gridxx  also  does  not  consider  the  effect  of  manmade  or  natural  levees 
nor  the  actual  flow  of  water  across  the  complex  landscape. 

Data  are  also  provided  for  the  results  of  the  cost  distance  grid  using  the 
previous  intermediate  flood  extent  (Figure  18).  This  grid  theme  shows  the 
shortest  path  distance  from  each  cell  back  to  the  closest  water  elevation  source 
point. 

In  Figure  19,  a  second  water  surface  grid  theme  is  generated  in  the  diagnostic 
mode  and  is  illustrated.  This  second  surface  is  the  result  of  all  the  gage, 
centerline  channel,  secondary  channel,  and  flow  path  water  surface  elevations. 
This  water  surface  grid  is  then  compared  to  the  DEM,  which  results  in  the  theme 
‘Improved  Surf  as  shown  in  Figure  20. 

Using  the  second  flood  event  approximation,  ‘Improved  Surf,  water 
connectivity  is  evaluated  and  a  final  flood  surface  is  generated  (Figure  21). 
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7  FEAT  Algorithms 


The  following  is  a  step-by-step  description  of  the  algorithms  and  processes 
used  to  generate  a  flood  surface.  Much  of  the  terminology  used  in  this  chapter  is 
compatible  with  the  ArcView  Avenue  programming  language. 

The  first  part  of  the  simulation  collects  and  evaluates  all  the  source  data  from 
the  user.  Typically  these  are  river  gage  water  elevations,  river  center  lines, 
secondary  channels,  and  river  channel  cross  sections.  Water  elevations  from  the 
gage  theme  are  associated  with  the  river  center-line  polyline  endpoints.  Once  the 
endpoints  have  elevation  values,  each  of  the  vertices  of  the  polyline  is  assigned 
an  elevation  based  on  the  beginning  and  endpoint  elevation  slope  and  distance 
from  the  beginning  of  the  polyline.  The  distance  calculated  is  the  distance  along 
the  polyline,  not  a  straight  line  distance.  The  vertices  attributed  to  the  elevation 
are  combined  with  the  gage  water  elevations  and  produce  a  temporary  point 
theme.  This  process  is  repeated  similarly  for  the  secondary  channels,  except  that 
the  elevation  for  all  the  vertices  is  assigned  by  finding  the  closest  center-line 
vertex  to  one  of  the  secondary  channel  endpoints.  The  elevation  from  this  vertex 
is  assigned  to  all  the  vertices  of  the  secondary  channel  and  the  vertices  are 
appended  to  the  temporary  point  theme. 

Once  the  point  theme  has  been  generated  with  all  the  necessary  water  surface 
elevations,  a  water  stage  surface  elevation  grid  is  generated.  This  water  surface 
grid  is  generated  at  a  large  cell  resolution  (100  m)  using  the  inverse  distance 
weighted  average  algorithm  using  a  minimum  of  12  sample  points.  A  large 
resolution  was  selected  at  this  time  to  speed  up  the  interpolation  process  and 
because  this  approximated  surface  only  provides  a  rough  estimate  of  the  extent  of 
the  flood  event. 

The  approximated  water  surface  is  now  compared  to  the  ground  elevation 
model  to  evaluate  a  first  approximated  flood  area.  If  the  elevation  of  the  water 
surface  is  greater  than  or  equal  to  the  ground  surface,  then  the  grid  cell  is 
flooded;  otherwise,  the  grid  is  not  flooded.  The  resulting  flood  area  does  not 
contain  water  surface  elevations.  It  only  indicates  all  grids  that  could  be  flooded 
under  these  stage  conditions.  It  does  not  consider  effects  of  man-made  or  natural 
levees  or  the  actual  flow  of  water  across  the  flood  plain. 

The  next  step  determines  the  effects  of  topography  on  overland  flow.  The 
initial  flood  surface  is  modified  to  eliminate  all  flooded  grid  cells  that  have  no 
direct  connectivity  to  any  of  the  input  stage  data  locations.  This  is  accomplished 
in  ArcView®  Spatial  Analyst  Extension  by  using  the  function  CostDistanceQ.  The 
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function  CostDistanceQ  calculates  the  least-  accumulative-cost  distance  over  a 
surface  to  a  set  of  source  points  and  also  produces  a  direction  surface  for  each 
cell  to  its  closest  source  point.  In  this  application  the  source  points  are  the  stage 
locations  in  the  main  river  channel,  vertices  of  the  center  lines,  and  vertices  of 
secondary  channels.  The  direction  surface  shows  the  overland  flow  path  from 
each  cell  to  the  closest  flood  water  source. 

Next  the  first  approximate  flood  surface  and  the  direction  surface  are  used  to 
determine  the  flood  surface  elevation  in  each  flooded  cell.  The  direction  surface 
is  sampled  at  user-defined  regular  intervals  and  analyzed  to  define  the  shortest 
water  flow  path  from  a  cell  back  to  the  source  point  across  the  flood  surface.  The 
water  surface  elevation  of  the  source  point  is  assigned  to  the  cell  and  all  cells  are 
traversed.  A  default  value  of  25  allows  a  25  by  25  sampling  grid  over  the  area. 
Since  this  is  a  uniform  sampling,  increasing  this  sampling  can  provide  for 
improved  results.  If  the  sampling  interval  is  too  large,  many  small  backwater 
features  (streams,  channels)  may  not  be  sampled  and  thereby  not  depict  flooding 
of  backwater  areas. 

Using  the  water  surface  elevation  of  the  sample  cells,  all  the  cells  in  the  flow 
path,  the  original  gage  elevations,  center  line,  and  secondary'  channel  vertices,  a 
second  approximate  flood  surface  is  interpolated.  This  surface  is  generated  at  the 
same  cell  size  as  the  ground  elevation  model  using  the  inverse  distance  weighted 
interpolation  method.  This  surface  is  generated  without  regard  to  topography. 

This  second  approximated  flood  surface  is  then  modified  to  eliminate  all 
flooded  grid  cells  that  do  not  have  direct  connectivity  to  any  of  the  input  water 
elevation  source  point  locations.  This  is  the  same  process  as  was  applied  to  the 
first  approximated  flood  surface  using  the  same  functions.  The  resulting  distance 
surface  is  converted  to  the  final  flood  area  surface. 

The  final  flood  area  surface  is  added  to  the  current  view  in  ArcView  and  the 
process  is  terminated. 
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Appendix  A 
Model  Listing 


DRIVER4.AVE 


Name:  LMK.Driver4 
Name:  Driver4 


Name:  Driver4.ave 


Purpose:  Test  Avenue  script  driver  for  user  input  of  line  and  point 
shape  files 


Calls:  LMK.line2dpts,  LM K.GetPointTheme,  LMK.Water-Flow 
LNK.streamline 

Comments:  This  is  different  from  Driver3  in  that  it  includes 
interpolated  water  surface  values  on  channel  centerlines 


Developer:  USA  ERDC-WES/  Mr.  Jerry  Ballard 
3909  Halls  Ferry  Road 
Vicksburg,  MS  39180 


Sponsor:  CELMK-ED-H  (Mr.  Ron  Goldman) 
4155  Clay  Street 
Vicksburg,  MS  39180-3435 


Notes: 
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'  3/25/1998  -  fixed  problem  with  deleting  the  direction  grid  from  first 
1  cost  distance  iteration.  For  some  reason  it  can  only  be  deleted 
'  before  it  is  used.  After  it  is  generated,  it  cannot  be  deleted 
'  from  the  same  script,  -jrb 

t 

'  3/25/1998  -  Added  the  item  name  to  the  final  flood  surface.  Not 

'  really  sure  what  will  happen  if  there  is  a  combined  point  and 

'  transect  inputs,  -jrb 
1 

'  10/15/1998  -  Rats!  ArcView3.1  upgrade  undid  my  fix  on  the  direction 

'grid. 

1 

'  10/15/1998  -  Added  the  LMK.steamline  to  interpolate  water  surface 
'  elevations  along  channel  centerlines  (just  like  FIST!). 

I 

'  10/21/1998  -  Added  the  LMK. Secondary  to  add  flat  water  surface 
'  elevations  along  secondary  channels. 

I 

theView  =  av.GetActiveDoc 


Get  a  list  of  Grid  Themes  and  have  user  select  one 


'  Get  a  list  of  all  the  Themes 

t 

IstAllThemes  =  theView.GetThemes 

IstGridThemes  =  list.Make 
1 

'  Get  a  list  of  all  available  Grid  Themes 

I 

for  each  thm  in  IstAllThemes 
if  ( thm.IS(GTheme))  then 
IstGridThemes.  Add(thm) 
end 
end 

'  Did  we  find  any? 

I 

if  ( IstGridThemes.Count  =  0 )  then 

MsgBox,Error("There  are  NO  Grid  Themes  in  this  View.", "Returning") 

retum(nil) 

end 

I 

'  Prompt  the  user  to  select  one  of  the  Grid  Themes 

I 

thmElevTheme  =  MsgBox.List(lstGridThemes, 

"Select  the  Grid  Theme  that  contains  the  terrain  elevations," 

"Terrain  Elevation  Grid") 
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Check  for  cancelation  of  procedure 


if  ( thmElevTheme  =  nil )  then 

retum(nil) 

end 


Get  a  list  of  Point  Themes  and  let  the  user  possibly 
select  one  for  processing. 


IstPointThemeInfo  =  av.Run("LMK.GetPointTheme",{}) 


If  a  point  theme  of  water  gages  was  selected,  then  we  can 
ask  for  the  theme  of  polylines  for  the  channel  center-line 
interpolation. 


if  ( IstPointThemeInfo.Count  >  1)  then 

IstInterpPointThemeInfo  =  av.Run("LMK.streamline",lstPointThemeInfo) 

'  thmlnterpPointTheme  =  lstlnterPointThemelnfo.Get(O) 

IstSecondPointThemeInfo  =  av.Run("LMK.  Secondary", IstInterpPointThemeInfo) 
end 


'  Get  a  list  of  Line  Themes  (cross-sectional  input)  and  let  the 
'  user  select  one  and  also  the  item  of  the  theme  to  use 


IstLineThemeInfo  =  av.Run("LMK.line2dpts",{}) 


'  Prompt  for  the  sampling  steps  for  use  in  LMK.Water-Flow 
'  Default  is  25 


strSample  Steps  = 

MsgBox.Input("  Specify  Sampling  Steps  for  Water  Flow  Calculation," 
"Specify  Sampling  Steps  for  Water  Flow  Calculation," 

"25") 


if  ( (strSampleSteps.IsNumber)  or  (strSampleSteps.AsNumber  >  0) )  then 

numSampleSteps  =  strSampleSteps.AsNumber 

else 

MsgBox.Error("Sampling  Steps  Provided  is  not  numeric.  Returning," 
"Error  in  Input") 
return  {nil} 
end 


'  Prompt  for  the  additive  constant  for  the  initial  watersurface 
'  Default  is  0 
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strAddElev  =  MsgBox,Input( 

"Specif)'  Additive  Constant  for  Initial  Surface  in  Meters," 
"Specify  Additive  Constant  for  Initial  Surface," 

"0") 


if  (  strAddElev.IsNumber )  then 
numAddElev  =  strAddElev.AsNumber 
else 

MsgBox.Error(" Additive  Constant  entered  is  not  numeric.  Returning," 
"Error  in  Input") 
return  {nil} 
end 


'  Some  convoluted  logic  here . 

f 

'  If  neither  a  point  nor  line  shape  file  was  selected,  then  cancel  the 
*  whole  process. 

I 

'  If  just  the  point  shapefile  was  selected,  then  continue  as  planned. 

'  If  just  the  line  shapefile  was  selected,  then  continue  as  planned. 

'  If  ♦BOTH*  were  selected,  then  we  have  to  merge  the  point  and  line 
'  shapefile  into  a  single  point  coverage.  UGH.  But  fortunately  the 
'  line  shapefile  was  just  converted  to  a  shapcpoint  file  so  we  can  just 
'  keep  adding  features  to  the  shapefile  and  then  continue. 


strFieldName  =  ""  'Kludge  to  keep  field  name  from  point  shapefile 

if  ( (IstPointThemeInfo. Count  <  3)  and  (IstLineThemeInfo.Count  <  2) )  then 
MsgBox.Infof'No  Theme  inputs  selected.  Ending  Procedure," 

"Procedure  Termination") 

return  nil 

end 

if  ( (IstPointThemeInfo. Count  >  1)  and  (IstLineThemeInfo.Count  <  2) ),  then 
Av.ShowMsg("Using  only  the  Point  Theme") 

I 

'  Copy/Clone  the  objects  into  required  variables 

t 

thmPointTheme  =  lstPoinfThemelnfo.Get(O) 
ftbPointTheme  =  thmPointTheme.GetFTab 
fldPointThemeStage  =  lstPointThemelnfo.Get(l) 
numPointThemeConv  =  lstPointThemeInfo.Get(2) 
fldPointThemeShape  =  ftbPointTheme.FindField("Shape") 
strFieldName  =  fldPointThemeStage.AsString 
end 
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if  ( (IstPointThemeInfo. Count  <  3)  and  (IstLineThemeInfo.Count  >  1) ),  then 

Av.ShowMsg(”Using  only  the  Line  Theme") 

1 

1  Copy  the  list  objects  into  required  variables 

I 

thmPointTheme  =  lstLineThemeInfo.Get(0) 
fldPointTheme  Stage  =  lstLineThemelnfo.Get(l) 
ftbPointTheme  =  thmPointTheme. GetFTab 
numPointThemeConv  =  1.0 

end 

if  (( IstPointThemeInfo. Count  >  1)  and  (IstLineThemeInfo.Count  >  1) ),  then 
Av.ShowMsg("Merging  line  and  point  themes") 

thmUserPointTheme  =  IstPointThemeInfo. Get(0) 
ftbUserPointTheme  =  thmUserPointTheme. GetFTab 
fldUserPointThemeStage  =  lstPointThemelnfo.Get(l) 
numUserPointConv  =  lstPointThemeInfo.Get(2) 

AdUserPointThemeShape  =  ftbUserPointTheme  .FindField("Shape") 

'  remember  actually  this  theme  is  a  point  theme  generated  from  the 
'  selected  line  theme 

thmUserLineTheme  =  lstLineThemelnfo.Get(O) 
fldUserLineThemeStage  =  lstLineThemelnfo.Get(l) 

ftbUserLineTheme  =  thmUserLineTheme.GetFTab 
fldUserLineThemeShape  =  ftbUserLineTheme.FindField("Shape") 

'  hopefully  we  have  all  the  info  now  to  merge  this  stuff 

I 

'  Read  though  the  point  shapefile  and  copy  the  records  into 
'  the  output  file. 

for  each  rec  in  ftbUserPointTheme 
shpUserPointObject  =  ftbUserPointTheme.RetumValue 
( AdUserPointThemeShape,  rec) 

numUserPointStage  =  ftbUserPointTheme. RetumValueNumber 
( fldUserPointThemeStage,  rec) 

numUserPointStage  =  numUserPointStage  *  numUserPointConv 
numPtRec  =  ftbUserLineTheme.  AddRecord 

ftbUserLineTheme. SetValueNumber(  fldUserLineThemeStage,  numPtRec, 
numUserPointStage) 

ftbUserLineTheme.  SetValue(  fldUserLineThemeShape,  numPtRec, 

shpUserPointObj  ect) 

end 
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'  Copy/Clone  the  objects  into  required  variables 

I 

thmPointTheme  =  thmUserLineTheme. Clone 
ftbPointTheme  =  ftbtiserLineTheme 
fldPointThemeStage  =  fldUserLineThemeStage 
fldPointThemeShape  =  AdUserLineThemeShape 
numPointThemeConv  =1.0 

end  '  End  of  conditional  using  both  line  and  point  themes 

I 

I********************************************************* 

I  %****%**%:*%%%%******* **%********* ************************ 

'  Additional  logic  for  adding/merging  line  and  point  themes 
'  specifically  for  the  center  line  interpolation 


'  (1)  if  a  line  theme  was  not  selected,  then  bail. 

if  ( lstInterpPointThemeInfo.Get(0).Is(  Theme  ) ),  then 

'  (2)  From  the  returned  list,  get  Theme,  Ftab,  and  Field 

thmlnterpPointTheme  =  IstlnterpPointThemelnfo.Get(O) 
ftblnterpPointTheme  =  thmlnterpPointTheme.GetFTab 
ftblnterpPoinflheme.Flush 

fldlnterpPoinflhemeStage  =  ftbInterpPointTheme.FindField("Stage") 
fldlnterpPointThemeShape  =  ftbInterpPointTheme.FindField("  Shape") 

'  (3)  The  previous  section  of  code  provides  a  point  theme 
'  ready-to-go.  We  really  aren't  sure  what  this  point 
'  theme  is,  so  we  need  to  append  it  to  ours. 

'  Important  variables: 

'  thmPointTheme  -  Point  Theme 
1  ftbPointTheme  -  FTab  of  point  theme 
'  fldPointThemeStage  -  Field  of  data  values 

'  fldPointThemeShape  -  Shape  of  the  Theme  (it  really  is  important) 

» 

'  Step  through  each  record  of  the  previous  point  theme 
for  each  rec  in  ftbPointTheme 
'  Get  the  shape 

shpPointTheme  =  ftbPointTheme  ,RetumValue(fidPointThemeShape, rec) 

'  Get  the  value  (number) 

numPointTheme  =  ftbPointTheme.RetumValueNumber(fldPointThemeStage,rec) 
numPointTheme  =numPointTheme  *  numPointThemeConv 

'  Append  a  record  to  the  interpolated  point  file 
numPtRec  =  ftblnterpPointTheme.AddRecord 
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'  Set  the  shape 

ftblnterpPointTheme .  SetValue(fldInterpPointThemeShape,  numPtRec, 
shpPointTheme) 

'  Set  the  value 

ftblnterpPointTheme. SetValueNumber(fldInterpPointThemeStage,  numPtRec, 

numPointTheme) 

end 


'  Process  the  secondary  channels  if  any 

t 

if  ( lstSecondPointThemeInfo.Get(0).Is(  Theme ) )  then 
thmSecondPointTheme  =  lstSecondPointThemeInfo.Get(0) 
ftbSecondPointTheme  =  thmSecondPointTheme. GetFTab 
ftbSecondPointTheme  .Flush 

AdSecondPointThemeStage  =  ftbSecondPointTheme.FindField("Stage") 
AdSecondPointTheme  Shape  =  ftbSecondPointTheme.FindField("Shape") 
for  each  rec  in  ftbSecondPointTheme 
'  Get  the  shape 
shpSecondPointTheme  = 

ftbSecondPointTheme.RetumValue(fldSecondPointThemeShape,rec) 

'  Get  the  value  (number) 
numSecondPointTheme  = 

ftbSecondPointTheme.RetumValueNumber(fIdSecondPointThemeStage,rec) 

'  numSecondPointTheme  =  numSecondPointTheme  *  numPointThemeConv 

'  Append  a  record  to  the  interpolated  point  file 
numPtRec  =  ftblnterpPointTheme.AddRecord 

'  Set  the  shape 

ftblnterpPointTheme.  SetValue(fldInterpPointThemeShape,  numPtRec, 
shpSecondPointTheme) 

'  Set  the  value 

ftblnterpPointTheme.  SetValueNumber(fldInterpPointThemeStage,  numPtRec, 
numSecondPointTheme) 

end 

end 

'  (4)  That  should  do  it.  Now  we  need  to  copy/clone  the 
'  objects  for  consistency. 

I 

'  Copy/Clone  the  objects  into  required  variables 

I 

thmPointTheme  =  thmlnterpPointTheme. Clone 
ftbPointTheme  =  ftblnterpPointTheme 
fldPointTheme  Stage  =  fldlnterpPointThemeStage 
fldPointThemeShape  =  fldlnterpPointThemeShape 
numPointThemeConv  =  1.0 
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end  '  of  center  of  streamline  interpolation  points 


'  Ok,  We  managed  to  navigate  though  the  mess  of  getting  or  merging  all 
'  necessary  shapefiles.  Important  variables  up  to  this  point: 

'  *  IMPORTANT  VARIABLES* 

I 

*  thmPointTheme  -  the  point  shape  file  for  determine  water  surface  elevation 
'  ftbPointTheme  -  the  FTab  of  that  shapefile 

*  fldPointThemeStage  -  the  item  name  or  field  containing  the  numeric  values 
'  of  the  water  elevation 

'  numPointThemeConv  -  conversion  factor  for  the  units  of  water  elevation 

I 

'  thmElevTheme  -  the  Grid  (GTheme)  Theme  containing  the  ground  (terrain) 
'  elevation  values  in  meters 


'  Make  ArcView  dump  everything  in  memory  about  ftbPointTheme  to  disk 
'  this  may  or  may  not  be  necessary'  depending  on  how  we  got  this  point  shapefile 

ftbPointTheme  .Flush 

Av.ClearMsg 


1  Prepare  for  the  initial  interpolation 


thePij  =  theView.GetProjection 


'  Access  the  terrain  elevation  grid  theme 


grdTopoElev  =  thmElevTheme  .GetGrid 


'  Do  all  initial  interpolation  at  a  large  cell  size  to  save  time 
1  (we'll  use  the  correct  cell  size  at  the  final  interpolation) 
numCellSize  =  100 

'numCellSize  =  grdTopoElev.GetCellSize 
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Set  analysis  area 


Grid.SetAnalysisCellSize(#GRID_ENVTYPE_VALUE,numCellSize) 
Grid.SetAnalysisExtent(#GRID_ENVTYPE_VALUE,  theView.RetumExtent) 


Convert  the  Water  Surface  Point  theme  to  a  grid 

This  converted  Grid  theme  will  be  used  later  for  calculating  cost  and 
flow  direction 


fldlDPointTheme  =  ftbPointTheme.FindField("Id") 
grdPointTheme  =  Grid.MakeFromFTab(ftbPointTheme,  thePq, 
fldlDPointTheme,  NIL) 

gthmPointTheme  =  GTheme.Make(grdPointTheme) 
gthmPointTheme.SetName("Gauge  Grid") 


I 

'  Set  the  radius  and  interpolation  for  surface  generation 

I 


theRadius  =  Radius.MakeVariable(12,nil) 
thelnterp  =  Interp.MakeIDW(2,theRadius,nil) 


This  provides  a  water  surface  elevation  determined  by  least  linear  distance 


Interpolate  the  initial  surface 


grdlnitial  Surface  =  Grid.MakeByInterpolation( 
ftbPointTheme,  thePij,  fldPointThemeStage,  thelnterp,  NIL) 

if  ( grdlnitial  Surface  .HasError)  then 
MsgBox.Error(" Initial  Interpolation  Failed.", "Error") 
return  nil 
end 


'  Access  the  terrain  elevation  grid  theme 
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grdTopoElev  =  thmElevTheme.GetGrid 


'  Determine  the  initial  flooded  surface 


'  Note:  Don't  forget  to  convert  vertical  units  here! ! 


grdRoughlnitFlood  = 

grdTopoElev  <=  ((grdlnitialSurface  *  numPointThemeConv.AsGrid)  + 

numAddElev .  AsGrid) 

if  ( grdRoughlnitFlood.HasError )  then 

MsgBox.Error("Initial  Flood  Delination  Failed.", "Error") 

return  nil 

end 


'  Get  rid  of  the  stupid  zeros  (sigh) 


thmRoughlnitFlood  =  GTheme.Make(grdRoughInitFlood) 
grdQGrid  =  grdRoughlnitFlood.  ExtractByAttributes("([value]  =  1)") 

thmlnitSurface  =  GTheme.Make(grdInitialSurface) 
thmQGrid  =  GTheme.Make(grdQGrid) 

thmlnitSurface. SetName("InitSurf') 

thmlnitSurface  .SetComments(fldPointThemeStage  .AsString) 

theView.AddTheme(thmlnitSurface) 

theView.AddTheme(thmQGrid) 

the  V  iew.  AddTheme  (thmPointTheme) 


'  Determine  all  water  areas  connected  to  main  line  stream 
'  -  this  effectively  gets  rid  of  floods  in  protected  areas 

t 

pq'Project  =  av.GetProject 

flnWorkDir  =  prjProject.GetWorkDir 

flnDirGrid  =  FileName.Merge(flnWorkDir.AsString,"dirgrid") 
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if  ( File.Exists(flnDirGrid) )  then 

if  ( Grid.DeleteDataSet(flnDirGrid) )  then 

'MsgBox.Info("Grid  Deleted! 

else 

'MsgBox.Info("Grid  *not*  deleted.","") 

end 

end 


Calculate  an  improved  water  surface  by  calculating  water  surface  elevation 
by  least  path  distance  (not  necessarily  a  linear  distance) 


Don't  be  confused  here:  I!!!!!!!!!!!!!!!!!!!!!!! 

grdQGrid  -  has  1  where  the  water  is  flooded.  This  is 
needed  for  CostDistance() 

grdlnitialSurface  -  has  the  water  surface  elevations  for 
tire  entire  area  extent. 


Calculate  cost  surface  and  direction  grid 


grdCostSurface  =  grdPointTheme.CostDistance(grdQGrid,flnDirGrid,NIL,NIL) 

thmCostSurface  =  GTheme.Make(grdCostSurface) 
thrnCostSurface  .SetName("Cost  Surface") 

theView.AddTheme(thmCostSurface) 


Using  the  direction  grid,  calculate  water  flow  paths 

Call  LMK.Water-Flow  with  Direction  Grid,  Initial  Surface,  ftbPointTheme, 
and  fldPointThemeStage 

This  returns  an  FTheme  with  point  water  surface  elevations. 
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thmSurfPts  =  av.Run("LMK.Water-Flow",  {flnDirG  rid, grdlniti  a]  Surface. 
ftbPointTheme.  fldPointThemeStage, 
numSampIeSteps} ) 


'  Try  to  delete  the  nasty  DirGrid 

I 

'if  ( File.Exists(flnDirGrid) )  then 

'  if  ( File.CanDelete(flnDirGrid) )  then 
'  Grid.DeleteDataSet(flnDirGrid) 

'  else 

'  MsgBox.Errorf'Cannot  remove  temporary  grid  file  "+flnDirGrid.AsString, 
'  "Unable  to  continue") 

'  return  nil 
'  end 
'end 


'  See  if  the  Water-Flow  script  did  its  job 
» 

if  ( thmSurfPts  =  Nil )  then 

MsgBox.Error("Flow  Direction  Calculation  Failed.", "Stopping  Procedure") 

exit 

end 


ftbSurfPts  =  thmSurfPts.GetFTab 
ftbSurfPts.Flush 

'theView.AddTheme(thmSurfPts) 

thmSurfPts.Invalidate(True) 

Grid.SetAnalysisCellSize(#GRID_ENVTYPE_VALUE,numCellSize) 
Grid.SetAnalysisExtent(#GRID_ENVTYPE_VALUE,  theView.RetumExtent) 


'  Interpolate  new  surface  with  the  surface  points 

♦ 

fldSurfPtsStage  =  ftbSurfPts.FindFieId("Stage") 


'  Set  the  radius  and  interpolation  for  surface  generation 


theRadius  =  Radius.MakeVariabIe(12,nil) 
thelnterp  =  Interp.MakeIDW(2,theRadius,nil) 
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Interpolate  the  second  water  surface  surface 


grdSecondSurface  =  Grid.MakeByInterpolation( 
ftbSurfPts,  thePij,  fldSurfPtsStage,  thelnterp,  NIL) 

if  ( grdSecondSurface. HasError)  then 
MsgBox.Error(" Second  Interpolation  Failed.", "Error") 
return  nil 
end 


'  Determine  new  flooded  surface 

t 

'  Change  cell  size  from  50  to  elevation  grid  cell  size 
1 

numCellSize  =  grdTopoElev.GetCellSize 


'  Set  analysis  area 
1 

Grid.SetAnalysisCellSize(#GRID_ENVTYPE_VAHJE,numCellSize) 
Grid.SetAnalysisExtent(#GRID_ENVTYPE_VALUE,  theView.RetumExtent) 

» 

1  Note:  Don't  forget  that  the  elevation  file  vertical  units  are  in  feet! 

grdSecondFlood  =  grdTopoElev  <=  (grdSecondSurface  * 

numPointThemeConv.AsGrid) 

if  ( grdSecondFlood.HasError )  then 

MsgBox.Error(" Second  Flood  Delination  Failed.", "Error") 

return  nil 

end 


'  Get  rid  of  the  stupid  zeros  (sigh) 


thmSecondFlood  =  GTheme.Make(grdSecondFlood) 

grdCleanSecondFlood  =  grdSecondFlood.ExtractByAttributes("([value]  =  1)") 

thmCleanSecondFlood  =  GTheme.Make(grdCleanSecondFlood) 
thmCleanSecondFlood.SetName("Flooded  Area") 

thmSecondSurface  =  GTheme.Make(grdSecondSurface) 
thmSecondSurface.SetName("Water  Surface") 
theView.AddTheme(thmSecondSurface) 
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thmSecondFlood.SetName("Improved  Surf') 

theView,AddTheme(thmSecondFlood) 

theView.AddTheme(thmCleanSecondFlood) 


'  Calculate  the  second  improved  water  surface,  by  calculating  water  surface 
elevation 

'  by  least  path  distance  (not  necessarily  a  linear  distance) 


'  Calculate  cost  surface  only  (We  don't  need  the  direction  grid  anymore) 

I 

flnWorkDir  =  pijProject.GefWorkDir 

flnNewDirGrid  =  FileName.Merge(flnWorkDir.AsString,"dirgridl") 

'if  ( File.Exists(flnNewDirGrid) )  then 
'  if  ( Grid.DeleteDataSet(flnNewDiiGrid) )  then 
'  MsgBox.Info("Grid  Deleted!!!","") 

'  else 

'  MsgBox.lnfo("Grid  *not*  deleted.","") 

'  end 
'end 


grdCostSurface2  = 

grdPointTheme.CostDistance(grdCleanSecondFlood,NIL,NIL,NIL) 

thmCostSurface2  =  GTheme.Make(grdCostSurface2) 
thmCostSurface2.SetName("Cost  Surface2") 
thmCostSurface2.SetActive(true) 

■the  View.  AddTheme(thmCostSurface2) 


Av.ShowMsg("Extracting  by  Attributes") 
grdFinalFlood  =  grdCostSurface2  >  O.AsGrid 
thmFinalFlood  =  GTheme.Make(grdFinalFlood) 
if  ( strFieldName.IsNull )  then 

thmFinalFlood. SetName("Flood  Surface:  H+fldPointThemeStage.AsString) 
else 

thmFinalFlood.SetName("Flood  Surface:  "+strFieldName) 
end 
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theView.AddTheme(thmFinalFlood) 


LMK.GetPointTheme.AVE 

'  Name:  LMK.GetPointTheme 
theView  =  Av.GetActiveDoc 


'  Get  a  list  of  all  Themes 

I 

IstAllThemes  =  theView.GetThemes 


'  Make  a  list  to  hold  all  feature  themes 
» 

IstFeatureList  =  list.Make 


'  Find  the  point  FThemes 
! 

for  each  th  in  IstAllThemes 

if  ( th.IS(FTheme))  then 

ftbUTheme  =  th.GetFTab 

shapefield  =  ftbUTheme.FindField("  Shape") 

shapetype  =  shapefield.GetType 

if  ( shapetype  =  #FIELD_SHAPEPOINT )  then 

lstFeatureList.Add(th) 

end 

end 

end 


'  Check  to  make  sure  the  list  is  not  empty 
» 

if  ( IstFeatureList.  Count  =  0 )  then 

MsgBox.Error("There  are  NO  Point  Feature  Themes  in  this  View'V'Error") 
return  {nil} 

end 

thmPointTheme  =  MsgBox.List(lstFeatureList, 

"Select  the  Point  Feature  Theme  that  contains  river  surface  elevations," 
"River  Surface  Point  Theme") 

if  ( thmPointTheme  =  nil ),  then 
return  {nil} 
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end 


'  Prompt  for  the  field  that  contains  water  elevations 

I 

ftbPointTheme  =  thmPomtTheme.GetFTab 
IstAllPointFields  =  ftbPointTheme.GetFields 
IstPointFields  =  List.Make 

for  each  fid  in  IstAllPointFields 

if  ( fld.IsTypeNumber ),  then 

lstPointFields.Add(fld) 

end 

end 

if  ( IstPointFields.Count  =  0  ),  then 

MsgBox.Error("There  are  NO  numeric  attributes  attached  for  the  points. ".''Error. 

Returning") 

return  {nil} 

end 

fldPointStage  =  MsgBox.List(lstPointFields, 

"Select  field  that  contains  the  river  elevations," 

"Select  Field") 


'  Make  sure  a  field  was  selected 

I 

if  ( fldPointStage  =  nil ),  then 

MsgBox.Error("No  Field  Selected.  Returning.", "Error  in  Field  Selection") 

return  {nil} 

end 


'  Make  sure  the  field  is  numeric 

I 

if  ( fldPointStage.IsTypeNumber.Not ),  then 

MsgBox.Error("Field  selected  is  not  a  numeric  field.  Ending  Procedure," 
"Error:  Non-numeric  field  for  point  theme") 
re  turn  {nil} 
end 


'  Ask  for  the  units 

boolUnits  =  MsgBox.YesNo("Are  the  elevations  in  feet?","Elevation  Units",true) 


'  Set  the  conversion  factor 
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if  ( boolUnits  =  true  ),  then 
numConvPointTheme  =  0.3048 
else 

numConvPointTheme  =  1.0 
end 

return  {  thmPointTheme,  fldPointStage,  numConvPointTheme  } 


LMK.Streamline.AVE 

'  Name:  LMK. streamline 
'  Name:  LMK.streamline 

'  Purpose:  Convert  line  shapelfile  into  a  densified  point  shape  file  with 
'  interpolated  water  elevation  attributes 

'  Author:  Jerry  Ballard 
'  USACE  Waterways  Experiment  Station 
'  Vicksburg,  MS  39180 

'  Sponsor:  CELMK-ED-H  (Mr.  Ron  Goldman) 

'4155  Clay  Street 
'  Vicksburg,  MS  39180-3435 

'  Initial  Coding:  October  14,1998 


Arguments:  (a  list) 

0  thmPointTheme  -  Gage  Point  Theme 

1  fldPointTheme  -  Water  Elevation  Field  to  use  from  Gage  Point  Theme 

2  numPointThemeConv  -  Conversion  Factor  (feet  to  meters,  m  to  ft.) 

Returns:  (a  list) 

0  thmlnterpPointTheme  -  Interpolated  Gage  Point  Theme 
Calls: 

Is  Called  By:  LMK.DriverX.ave 


'  Get  current  view 

I 

theProject  =  Av.GetProject 
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theView  =  Av.GetActiveDoc 


if  ( theView  =  nil )  then 

MsgBox.Info("The  View  was  not  found.  Retuming","Error") 

return  {nil} 

end 


'  Get  point  field  info  from  past  arguments 


thmPointTheme  =  Self.Get(O) 
ftbPointTheme  =  thmPointTheme.GetFTab 
fldPointThemeStage  =  Self.Get(l) 
numPointThemeConv  =  Self.Get(2) 


'  Get  a  list  of  all  Themes 
1 

IstAllThemes  =  theView.GetThemes 


'  Make  a  list  to  hold  all  feature  Themes 

I 

IstFeatureList  =  list.Make 


'  Find  the  line  FThemes 

for  each  th  in  IstAllThemes 

if  ( th.IS(FTheme))  then 

ftbUTheme  =  th.GetFTab 

shapefield  =  ftbUTheme.FindField("Shape") 

shapctype  =  shapefield.GetType 

if  ( shapetype  =  #FIELD_SHAPELINE )  then 

IstFeatureList.  Add(th) 

end 

end 

end 


'  Check  to  make  sure  the  list  is  not  empty 

I 

if  ( IstFeatureList. Count  =  0  )  then 

MsgBox.Error("There  are  NO  Line  Feature  Themes  in  this  View'V'Error") 

retum{nil} 

end 

thmLineTheme  =  MsgBox.List(lstFeatureList 

"Select  the  Line  Feature  Theme  that  contains  channel  centerlines," 
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"Channel  Surface  Line  Theme") 

if  ( thmLineTheme  =  nil ),  then 

return  {nil} 

end 


'  Prompt  for  the  field  that  contains  starting  gage  id 

I 

ftbLineTheme  =  thmLineTheme.GetFTab 
IstLineFields  =  ftbLineTheme  .GetFields 

if  ( IstLineFields. Count  <  3  ),  then 

MsgBox.Error("There  are  NO  attributes  attached  for  the  lines,"  "Error, 

Returning") 

retum{nil} 

end 

fldStartLineStage  =  MsgBox.List(lstLineFields, 

"Select  field  that  contains  the  starting  gage  ids," 

"Select  Field") 


'  Make  sure  a  field  was  selected 

t 

if  ( fldStartLineStage  =  nil ),  then 

MsgBox.Error("No  Field  Selected.  Returning. ","Error  in  Field  Selection") 

retum{nil} 

end 

fldEndLine  Stage  =  MsgBox.List(lstLineFields, 

"Select  field  that  contains  the  ending  gage  ids," 

"Select  Field") 


'  Make  sure  a  field  was  selected 
if  ( fldEndLineStage  =  nil ),  then 

MsgBox.Error("No  Field  Selected.  Returning.  ","Error  in  Field  Selection") 

return  {nil} 

end 

strDensity Interval  =  MsgBox.Input("  Specify  Density  Interval  for  Stream 
Centerlines  (meters)," 

"Density  Interval","  1000") 

if  ( strDensity  hitcrval.IsNumber ),  then 
numDensity  Interval  =  strDensitylntervaLAsNumber 
else 
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MsgBox,Error("Interval  Specified  is  not  numeric.  Retuming.","Error  in  Density 

Interval") 

return  {nil} 

end 


'  Specify  the  output  shapefile... 


'  Create  a  temporary  shape  file  for  the  water  points 

I 

flnCWD  =  theProject.GetWorkDir 
shpName  =  flnCWD.MakeTmp("cline","shp") 

'defaultName  =  FileName.Make("$HOME").MakeTmp(''shape","shp'') 

'shpName  =  FileDialog.Put(  defaultName,"*. shp","Output  Shape  File"  ) 

if  (shpName  =  nil)  then 

exit 

end 

shpName .  SetExtension("  shp ") 

ftbPointFTab  =  Ftab.MakeNew(shpName,Point) 

IstFields  =  List.Make 

lstFields.Add(Field.Make("Stage",#FIELD_FLOAT,  16,2)) 
ftbPointFTab  .AddFields(lstFields) 
fldPointShape  =  ftbPointFTab.  FindField("Shape") 
fldPointStage  =  ftbPointFTab.  FindField("  Stage") 


'  Process  each  line 


fldLineShape  =  ftbLineTheme.FindField("Shape") 
fldPointThemelD  =  ftbPointTheme.FindFieId("Id") 

for  each  rec  in  ftbLineTheme 

1  Get  Starting  and  Ending  Stage  IDs  for  lookup 

f 

numStartLineStage  =  ftbLineTheme.RetumValueNumber(  fldStartLineStage, 
rec) 

numEndLineStage  =  ftbLineTheme.RetumValueNumber(  fldEndLineStage,  rec) 


'  Search  the  point  gage  FTab  and  find  the  actual  water  level  at  each  stage 

I 

for  each  srec  in  ftbPointTheme 

numID  =  ftbPointTheme.RetumValueNumber(  fldPointThemelD,  srec ) 
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if  ( numID  =  numStartLineStage  ),  then 
numStartWaterStage  =  ftbPointTheme.RetumValueNumber( 
fldPointThemeStage,  srec ) 

numStartWaterStage  =  numStartWaterStage  *  numPointThemeConv 
end 

if  (  numID  =  numEndLine  Stage  ),  then 
numEndWaterStage  =  ftbPointTheme.RetumValueNumber( 
fldPointThemeStage,  srec ) 

numEndWaterStage  =  numEndWaterStage  *  numPointThemeConv 
end 

end 

t 

'  Get  the  shape  object  (it  should  be  a  Poly  Line) 

! 

shpLineObject  =  ftbLineTheme.RetumValue(  fldLineShape,  rec ) 

I 

'  Density  at  an  interval 

I 

plineDensified  =  shpLineObject. RetumDensified(numDensitylnterval) 

» 

'  Get  the  length 
! 

numPlineLength  =  plineDensified.RetumLength 
» 

'  Compute  water  slope 

t 

numWaterSlope  =  ( numEndWaterStage  -  numStartWaterStage  )  / 
numPlineLength 

'msgbox.report(numPlineLength.AsString  +  " "  +  numWaterSlope.AsString,"") 


'  Convert  the  polyline  to  a  list  to  extract  points 
» 

IstPolyLine  =  plineDensified.AsList 

for  each  shapePart  in  IstPolyLine 
for  each  xyPoint  in  shapePart 

numPosition  =  plineDensified.PointPosition(  xyPoint )  / 100.0 
numWaterLevel  =  numStartWaterStage  +  (numWaterSlope  *  numPosition  * 
numPlineLength ) 

numPtrec  =  ftbPointFTab.AddRecord 

ftbPointFTab.SetValueNumber(  fldPointStage,  numPtrec,  numWaterLevel ) 

ftbPointFTab.SetValue(  fldPointShape,  numPtrec,  xyPoint ) 

end 

end 

end 
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thmlnterpPointTheme  =  FTheme,Make(ftbPointFTab) 
return  {thmlnterpPointTheme} 


LMK.Secondary.AVE 

'  Name:  LMK. Secondary 

I 

'  Arguments:  (a  list) 

’  0  thmPointTheme  -  Gage  Point  Theme 

'  1  fldPointTheme  -  Water  Elevation  Field  to  use  from  Gage  Point  Theme 
'  2  numPointThemeConv  -  Conversion  Factor  (feet  to  meters,  m  to  ft.) 

I 

'  Get  current  view 

I 

theProject  =  Av.GetProject 
theView  =  Av.GetActiveDoc 
if  ( theView  =  nil )  then 

MsgBox.Info("The  View  was  not  found.  Retuming'V'Error") 

return  {nil} 

end 


'  Get  point  field  info  from  past  arguments 

I 

thmPointTheme  =  Self.Get(O) 
ftbPointTheme  =  thmPointTheme. GetFTab 
fldPointThemeStage  =  ftbPointTheme.FindField("Stage") 


'  Get  a  list  of  all  Themes 
1 

IstAIIThemes  =  theView.GetThemes 


'  Make  a  list  to  hold  all  feature  Themes 

I 

IstFeatureList  =  list.Make 


'  Find  the  line  FThemes 

t 

for  each  th  in  IstAIIThemes 

if  ( th.IS(FTheme)),  then 

ftbUTheme  =  th.GetFTab 

shapefield  =  ftbUTheme.FindField("Shape") 

shapetype  =  shapefield.GetType 
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if  ( shapetype  =  #FIELD_SHAPELINE  )  then 

lstFeatureList.Add(th) 

end 

end 

end 


'  Check  to  make  sure  the  list  is  not  empty 
» 

if  ( IstFeatureList.Count  =  0  ),  then 

MsgBox.Error("There  are  NO  Line  Feature  Themes  in  this  View","Error") 

return  {nil} 

end 

thmLineTheme  =  MsgBox.List(lstFeatureList, 

''Select  the  Line  Feature  Theme  that  contains  secondary  channels," 
"Secondary  Channel  Line  Themes") 

if  ( thmLineTheme  =  nil ),  then 

return  {nil} 

end 


ftbLineTheme  =  thmLineTheme.GetFTab 

fldLineShape  =  ftbLineTheme  .FindField("Shape") 

fldPointShape  =  ftbPointTheme.FindField("Shape") 

» 

'  Create  a  temporary  shape  file  for  the  secondary  water  points 
\ 

flnWorkDir  =  theProject.GetWorkDir 
flnSecondPointTheme  =  flnWorkDir.MakeTmp("sline","shp") 

ftbSecondPointFTab  =  FTab.MakeNew(flnSecondPointTheme,Point) 
IstFields  =  List.Make 

lstFields.Add(Field.Make("Stage",#FIELD_FLOAT,  1 6,2)) 
ftbSecondPointFTab.AddFields(lstFields) 
fldSecondPointShape  =  ftbSecondPointFTab.FindField("Shape") 
fldSecondPointStage  =  ftbSecondPointFTab.  FindField("  Stage") 

'Process  each  line 

for  each  rec  in  ftbLineTheme 

shpLineObject  =  ftbLinetheme.RetumValue(  fldLineShape,  rec ) 

'MsgBox.Report(shpLineObject.AsString,"") 

plineDens  =  shpLineObject.RetumDensified(150) 

'IstPolyLine  =  shpLineObject.  AsList 
IstPolyLine  =  plineDens.AsList 
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for  each  shapePart  in  lstPolyLine 

xj'PointFirst  =  shapePart.Get(O) 

’MsgBox.Info(xyPointFirst.AsString,"") 

ftbPointTheme.SelectByPoint(xyPointFirst,100.0,#VTABjSELTYPE_NEW) 

if  ( ftbPointTheme.GetSelection.Count  >  0 )  then 
setPointTheme  =  ftbPointTheme.GetSelection 
for  each  prec  in  setPointTheme 

numPointStage  =  ftbPointTheme. Return  ValueNumber( 
fldPo  intTheme  Stage, prec) 

'MsgBox.Info(numPointStage.AsString,,M') 

end 

end 

for  each  xyPoint  in  shapePart 

numPtrec  =  ftbSecondPointFTab  .AddRecord 

ftbSecondPointFTab.SetValueNumber(  fldSecondPointStage,  numPtrec, 
numPointStage ) 

ftbSecondPointFTab. SetValue(  fldSecondPointShape,  numPtrec,  xyPoint ) 

'MsgBox.Info(xyPoint.AsString,"") 

end 

end 

end 

thmSecondPointTheme  =  FTheme.Make(ftbSecondPointFTab) 
ftbPointTheme  .Flush 

return  {thmSecondPointTheme,  1} 


LMK.line2dpts.ave 

'  Name:  LMK.line2dpts 
'  Name:  line2dpts.ave 

t 

'  Purpose:  Convert  line  shape  file  into  a  densified  point  shape  file  with 
'  water  elevation  attributes 


'  Get  current  project 

I 

theProject  =  Av.GetProject 
theView7  =  Av.GetActiveDoc 
if  ( theView  =  nil )  then 

MsgBox.Info("The  View  was  not  found.  Returning", "Error") 

return  {nil) 

end 
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Get  a  list  of  all  Themes 


IstAllThemes  =  theView.GetThemes 


'  Make  a  list  to  hold  all  feature  Themes 

I 

IstFeatureList  =  list.Make 


'  Find  the  line  FThemes 
» 

for  each  th  in  IstAllThemes 

if  ( th.IS(FTheme))  then 

ftbUTheme  =  th.GetFTab 

shapefield  =  ftbUTheme.FindField("Shape") 

shapetype  =  shapefield.GetType 

if  (  shapetype  =  #FIELD_SHAPELINE  )  then 

IstFeatureList .  Add(th) 

end 

end 

end 


'  Check  to  make  sure  the  list  is  not  empty 

l 

if  ( IstFeatureList.Count  =  0 )  then 

MsgBox.Error("There  are  NO  Line  Feature  Themes  in  this  View", ''Error") 

return  {nil} 

end 

thmLineTheme  =  MsgBox.List(lstFeatureList, 

"Select  the  Line  Feature  Theme  that  contains  river  cross  section  elevations", 
"River  Cross  Section  Line  Theme") 

if  ( thmLineTheme  =  nil )  then 

return  {nil} 

end 


'  Prompt  for  the  field  that  contains  water  elevations 
! 

ftbLineTheme  =  thmLineTheme.GetFTab 
IstLineFields  =  ftbLineTheme.GetFields 

if  ( lstLineFields.Count  <  3  )  then 

MsgBox.Error("There  are  NO  attributes  attached  for  the  lines. ","Error, 

Returning") 

return  {nil} 
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end 


fldLineStage  =  MsgBox.ListQstLmeFields, 
"Select  field  that  contains  the  river  elevations," 
"Select  Field") 


'  Make  sure  a  field  was  selected 

I 

if  ( fldLineStage  =  nil )  then 

MsgBox.Errorf'No  Field  Selected.  Returning. ","Error  in  Field  Selection") 

return  {nil} 

end 


'  Ask  for  the  units 

t 

boolUnits  =  MsgBox.YesNo("Are  the  elevations  in  feet?","Elevation  Units",true) 


'  Set  the  conversion  factor 

I 

if  ( boolUnits  =  true  )  then 

nuxnConvLineTheme  =  0.3048 
else 

numConvLineTheme  =  1.0 
end 


'  Prompt  for  new  shape  file  name 

I 

'  Specify  the  output  shapefile... 

! 

defaultName  =  FileName.Make("$HOME").MakeTmp("shape","shp") 

shpName  =  FileDialog.Put(  defaultName,"*  ,shp", "Output  Shape  File"  ) 

if  (shpName  =  nil),  then 

exit 

end 

shpName.SetExtension("shp") 

ftbPointFTab  =  Ftab.MakeNew(shpName, Point) 

IstFields  =  List.Make 

lstFields.Add(Field.Make("Stage",#FIELD_FLOAT,  16,2)) 
ftbPointFTab  AddFields(lstFields) 
fldPointShape  =  ftbPointFTab.FindField("Shape") 
fldPointStage  =  ftbPointFTab. FindField("  Stage") 


'  Process  each  line 
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fldLineShape  =  ftbLineTheme.FindField("  Shape") 
for  each  rec  in  ftbLineTheme 


'  Get  the  line  object  and  its  field  attribute 

I 

shpLineObject  =  ftbLineTheme.RetumValue(  fldLineShape,  rec ) 
numLineStage  =  ftbLineTheme. RetumValueNumber(  fldLineStage,  rec) 
numLineStage  =  numLineStage  *  numConvLineTheme 

'msgbox.report(shpLineObject.AsString,"") 

'  Density  the  line  to  my  liking  (20-m) 

t 

plineDensified  =  shpLineObject.RetumDensified(20) 

t 

'  Convert  the  polyline  to  a  list  to  extract  points 
! 

IstPolyLine  =  plineDensified.AsList 

for  each  shapePart  in  IstPolyLine 
for  each  xyPoint  in  shapePart 
numPtrec  =  ftbPointFTab.AddRecord 

ftbPointFTab.SetValueNumber(  fldPointStage,  numPtrec,  numLineStage  ) 

ftbPointFTab.SetValue(  fldPointShape,  numPtrec,  xyPoint ) 

end 

end 

'msgbox.reportfplineDensified.AsString,"") 

end 


thmPointTheme  =  FTheme.Make(ftbPointFTab) 
return  {thmPointTheme,  fldPointStage  } 

LMK.Water-Flow.ave 

'  Name:  LMK.Water-Flow 
» 

'  Name:  LMK.Water-Flow 

I 

'  Title:  Calculate  shortest  water  flow  path 

I 

'  Topics:  CostDistance,  CellValue 


Appendix  A  Model  Listing 


A27 


'  Description:  This  script  calculates  the  shortest  water  flow  path  of  a  flood 
'  surface  using  a  direction  grid  from  CostDistance. 

f 

'  Requires:  Spatial  Analyst 

f 

'  Self: 

'  (0)  File  name  of  the  direction  grid 
'  (1)  Grid  of  water  surface  elevations 
'  (2)  FTab  of  the  stage  point  theme 
'  (3)  Field  of  surface  elevations  of  the  stage  point  theme 
'  (4)  Number  of  Sample  Steps 

f 

'  Returns:  Theme  composed  on  points  with  attributed  elevations 
* 

'  Author:  Jerry  Ballard 
'  USACE-WES-EN-C 
1  3909  Halls  Ferry  Road 
'  Vicksburg,  MS  39180 

I 

'  Sponsor:  USACE  Mississippi  Valley  District 

I 

'  Date:  03  Nov,  1997 

l 

'  Modifications: 

t 

'  10/16/98  -  Added  number  of  sample  steps  to  be  passed  in  instead 
'  of  constant 


*  Check  to  make  sure  we  have  at  least  five  arguments 

f 

if  ((Self.IS(List).Not)  or  (Self  count  <  5))  then 

return  nil 

end 

flnDirGrid  =  Self.Get(O) 
grdhritialSurfa.ee  =  Self.Get(l) 
ftbPointTheme  =  Self.Get(2) 
fldPointTheme  =  Self.Get(3) 
numSampleSteps  =  Self,Get(4) 


‘thmSurfPts  =  av.Run("LMK.Water-Flow",(flnDirGrid,grdInitialSurface, 
*  ftbPointTheme,  fldPointTheme, 

'  numSampleSteps}  ) 
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'  Convert  direction  grid  filename  to  a  source  name 

srcDirGrid  =  Grid.MakeSrcName(flnDirGrid.AsString) 

'  Check  and  see  if  it  is  valid 
» 

if  (srcDirGrid  =  nil )  then 

MsgBox.Error("LMK.Water-flow  script  got  an  invalid  Direction  Grid 

Name","Error") 

return  nil 

end 


'  Convert  to  Theme  and  get  Grid  property 

f 

thmDirGrid  =  Theme.Make(srcDirGrid) 
grdDirGrid  =  thmDirGridGetGrid 


1  Create  a  temporary  shape  file  for  the  water  points 

I 

pijProject  =  av.GetProject 

flnWorkDir  =  pijProject  GetWorkDir 

'flnCWD  =  FileName  GetCWD 

flnWaterPoints  =  flnWorkDir.MakeTmp("wfp","shp") 

ftbWaterPoints  =  FTab.MakeNew(flnWatcrPoints, Point) 


'  Add  fields 
» 

IstFields  =  List.Make 

lstFields.Add(Field.Make("Stage",#FIELD_FLOAT,16,2)) 
ftbWaterPoints.AddFields(lstFields) 
fldShapeField  =  ftbWaterPoints. FindField("Shape") 
fldStageField  =  ftbWaterPoints.FindField("Stage") 


'  Copy  all  the  field  data  from  the  initial  stage  point  theme 
» 

AdPointThemeShapeField  =  ftbPointTheme.FindField("Shape") 
for  each  rec  in  ftbPointTheme 

shpCurrentShape  =  ftbPointTheme.RetumValue(  AdPointThemeShapeField,  rec) 
numCurrentStage  =  ftbPointTheme.  RetumValueNumber(  fldPointTheme,rec) 

xyPoint  =  shpCurrentShape. RetumCenter 
recWaterPoints  =  fitbWaterPoints.AddRecord 
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ftbWaterPoints.SetValucNumber(  fldStageField,  recWaterPomts, 
num  CurrentStage) 

ftbWaterPoints . SetVaJue(  fldShapeField,  recWaterPoints,  xyPoint) 
end 


'  Get  Grid  characteristics 

I 

numCellSize  =  grdlnitialSurface.GetCellSize 
rectGridExtent  =  grdlnitial  Surface  .GetExtent 
numRows  =  grdlnitialSurface.GetNumRowsandCols.Get(O) 
numCols  =  grdInitialSurface.GetNumRowsandCols.Get(l) 

thePij  =  Pij.MakeNull 

theX  =  rcctGridExtcntGetLeft 
theY  =  rectGridExtent.GetTop 


'  Begin  sampling  loop 

I 

numXinc  =  (rectGridExtent.  GetRight  -  rectGridExtent.GetLeft  ).abs  / 
numSampleSteps 

numYinc  =  (rectGridExtent.GetTop  -  rectGiidExtent.GetBottom  ).abs  / 
numSampleSteps 

av.ShowMsg("Sampling  water  flow  routes...") 

for  each  ystep  in  1  ..numSampleSteps 

doMore  =  av.SetStatus(  (ystep/numSampleSteps)  *  100 ) 

for  each  xstep  in  1.. numSampleSteps 

ptX  =  theX  +  (numXinc  *  (xstep  -  1))  +  ((numXinc  /  2.0)  *  ((ystep  -1)  mod  2)) 
ptY  =  theY  -  (numYinc  *  (ystep  -  1)) 

p  =  Point. Make(ptX.ptY) 

theCellValue  =  grdDirGrid.CelIValue(p,thePij) 

IstPointList  =  List.Make 
lstPointList.Add(p) 

for  each  n  in  1..500 

if  ( theCellValue  .IsNull),  then 
'  MsgBox.Info("Null  point  value",'"') 
break 
end 
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if  ( theCellValue  =  1 ),  then 
ptX  =  ptX  +  numCellSize 

elseif  ( theCellValue  -  2  ),  then 
ptX  =  ptX  +  numCellSize 
ptY  =  ptY  -  numCellSize 

elseif  ( theCellValue  =  3  ),  then 
ptY  =  ptY  -  numCellSize 

elseif  ( theCellValue  =  4  ),  then 
ptX  =  ptX  -  numCellSize 
ptY  =  ptY  -  numCellSize 

elseif  (theCellValue  =  5  ),  then 
ptX  =  ptX  -  numCellSize 

elseif  ( theCellValue  =  6 ),  then 
ptX  =  ptX  -  numCellSize 
ptY  =  ptY  +  numCellSize 

elseif  ( theCellValue  =  8  ),  then 
ptX  =  ptX  +  numCellSize 
ptY  =  ptY  +  numCellSize 

elseif  ( theCellValue  =  7),  then 
ptY  =  ptY  +  numCellSize 

elseif  ( theCellValue  =  0 ),  then 
'  Get  the  water  surface  elevation  value 
numWaterElev  =  grdInitialSurface.CellValue(p,thePij) 

1  write  out  all  the  points  to  the  shape  file 

for  each  mylstpt  in  IstPointList 

recWaterPoints  =  ftbWaterPoints.AddRecord 

ftbWaterPoints .  SetV  alueNumber( 

fldStageField,  recWaterPoints,  numWaterElev) 

ftbWaterPoints. SetValue(  fldShapeField,  recWaterPoints,  mylstpt) 

end 

1  Don't  forget  this! 

break 

end 

p  =  Point.Make(ptX,  ptY) 

theCellValue  =  grdDirGrid.CellValue(p,thePij) 

lstPointList.Add(p) 

end  '  each  n  loop 
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IstPointList.Empty 

end 

end 

av.ClearMsg 

thmWaterPoints  =  Ftheme.Make(ftbWaterPoints) 
return  thmWaterPoints 
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